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Abstract 

Entanglement entropy in the vacuum state of local field theories exhibits an area 
law. However, nonlocal theories at large N and strong coupling violate this area law. 
In these theories, the leading divergence in the entanglement entropy is extensive for 
regions smaller than the effective nonlocality scale and proportional to this effective 
nonlocality scale for regions larger than it. This raises the question: is a volume law a 
generic feature of nonlocal theories, or is it only present at strong coupling and large 
A? 

This paper investigates entanglement entropy of large regions in weakly coupled 
nonlocal theories, to leading order in the coupling. The two theories studied are </> 4 
theory on the noncommutative plane and q i > 4 theory with a dipole type nonlocal modi¬ 
fication using a fixed nonlocality scale. Both theories are found to follow an area law to 
first order in the coupling, hence no evidence is found for a volume law. This indicates 
that, perturbatively the nonlocal interactions considered are not generating sufficient 
entanglement at distances of the nonlocality scale to change the leading divergence, at 
least to first order in the coupling. An argument against volume laws at higher orders 
is also presented. 
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1 Introduction and Summary 

Entanglement entropy has recently attracted interest as a way to study the correlations 
between degrees of freedom in a quantum state. Local held theories generally exhibit what 
is know as an area law behaviour, where the leading divergence in the entanglement entropy 
of a spatial region is proportional to the area of the boundary of that region. That is, 
S ~ \dA\A d ~ 2 3 , where S is the entanglement entropy, \dA\ the area of the boundary of 
the region and A is the momentum scale of the UV regulator of the theory, for example the 
inverse of a lattice spacing]^] However, recent holographic studies of strongly coupled nonlocal 
theories have found a volume law behaviour instead MISM. That is, for a nonlocality 
scale Z, S ~ |v4| 1 for regions much smaller than Z and S ~ Z|<9v4|A 1 ^ -1 for regions much 
larger than Z [5]. Note that entanglement entropy of large regions is sufficient to differentiate 
this type of volume law from an area law, as the entanglement entropy is proportional to the 
length scale of the nonlocality times an additional factor of the UV regulator. To summarise, 

area law : S ~ \dA\A d ~ 2 , 

volume law : S ~ |A|A d_1 , (small regions) 

S ~ Z|cL4|A d_1 . (large regions) 

These results can be understood intuitively by assuming that all the degrees of freedom 
within the range of the nonlocality are equally entangled with each other. Then, for regions 
much smaller than Z, all the degrees of freedom inside the region, not only those near the 
boundary, are entangled with degrees of freedom outside. For regions much larger than Z, all 
the degrees of freedom within a distance Z of the boundary are entangled with those outside. 
In both cases, the number of degrees of freedom strongly entangled across the boundary is 
proportional to A d_1 rather than the A d ~ 2 expected from an area law. 

A natural question is whether this behaviour is generic to nonlocal theories or if it is con¬ 
fined to a strongly coupled, large N regime. One approach is to study entanglement entropy 
for a free scalar held on the fuzzy sphere 0 0191 HU] . This turns out to be proportional 
to the aresj^Jfor small polar caps [9J [TO] . However, two issues arise which question whether 
this should be characterised as a volume law. First, the dependence of the entanglement 
entropy on the UV regulator does not match the volume law described above. Second, the 
entanglement entropy does not scale like the number of degrees of freedom contained in the 
polar cap, as the degrees of freedom are not uniformly distributed across the sphere. Instead 
it scales as the number of degrees of freedom near the boundary 0i- Another limitation of 
this theory is that the nonlocality scale is tied to the size of the sphere so it is not possible 
to study regions much larger than the nonlocality scale. 

Another approach is to study a free held theory on a lattice with a nonlocal kinetic term, 
in which case a volume law was found mi- 

This paper investigates the role of interactions in this question by considering two theories 
with nonlocal interactions: scalar A</> 4 theory on the noncommutative plane and A0 4 theory 
with a dipole type nonlocal modification with fixed nonlocality scale. The leading divergence 
in entanglement entropy of large regions is calculated to leading order in perturbation theory 

2 See for example p] for a review of area laws in entanglement entropy. 

3 The fuzzy sphere is a 2 dimensional surface, thus |A| is an area and |cb4| is a circumference. 
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and is not found to be proportional to the length scale of the nonlocality, hence no evidence 
of a volume law is found. Instead, the leading divergence in both theories has the same 
form as the standard local Xcft 4 theory which follows an area law. This result indicates 
that, perturbatively these nonlocal interactions are not generating sufficient entanglement 
at distances of the nonlocality scale to change the leading divergence, at least to first order 
in the coupling. 

The free theory with A = 0 for both of these nonlocal theories is equivalent to the regular 
commutative A (f) 4 theory. There is no modification of the entanglement entropy at this order. 
Perturbation theory can be used to study the nonlocal theories at small A. 

The entanglement entropy is calculated using the replica trick and the formula S = 
~d n [1 n Z n — nln Zi] n=1 , where Z n is the partition function of the held theory defined on an 
n-sheeted space H 21 na mi- This partition function can be reduced to computing vacuum 
bubble diagrams and the 0(A) contribution in perturbation theory comes from bubble dia¬ 
grams with one vertex and two loops. Consistent with the results of previous investigations 
of perturbative noncommutative theories [15], the planar diagrams in the nonlocal theories 
give the standard commutative result, which is S G\ (0) f dxd n=1 G n (x ) ~ A1_lA 2 ln(A/m), 
where A± is the (infinite) area of the boundary of our region, A our UV regulator, m our 
IR regulator and G n is the Green’s function on the n-sheeted space used in the replica trick 
This contribution follows an area law, as S' oc Aj_A 2 up to logarithmic corrections. 

The nonlocality only affects the nonplanar diagram. This diagram contributes a term 
of the form S ~ Gi(0, Ax) f dxd n= iG n (x, x + Ax) ~ (Axp ffi/(A, m, Ax), where now Ax 
corresponds to a translation from the nonlocality. 

In the dipole theory, Ax is proportional to the fixed dipole length. Thus the nonplanar 
diagram has only a logarithmic IR divergence and is subleading compared to the planar 
diagrams. In the noncommutative theory the translation along the noncommuative plane is 
proportional to the momentum in the other noncommutative direction, so this contribution 
must be integrated over this momentum. If we don’t impose an IR regulator, the momentum 
controlling the translation is allowed to vanish and G(0, Ax) —> Gi(0) ~ A 2 . This gives a 
contribution that is of the same order as the planar diagrams. However, if we impose an 
IR regulator, Ax has a minimal value and this divergence can be reinterpreted as an IR 
divergence. This is familiar from the UV/IR connection described for example in [15] . 

Our results for the 0(A) contribution to the entanglement entropy, Si, are 

real scalar . SI SAiSpianar T AiSnonpianar 

charged scalar : S± =(2A 0 + XGS v u^ r + Ai Anon planar - 

where Spianar and SAnp ianar denote the contributions from planar and nonplanar diagrams 
respectively. 


2 



The leading divergences from these diagrams in each of the theories considered are 


Commutative theory : 

Noncommutative plane : 
Dipole theory : 


S', 


planar 


A ± A 2 

2 10 3 2 7t 3 


5, 


nonplanar ^planar 


= Sn 


In 


A 2 
4 m 2 


5, 


nonplanar 


A_lA 2 - ln ( 
2 9 3 2 7t 3 1 - 


^nonplanar IS Subleadillg, 


9 2 m 2 A 2 ^ 

9 2 m 2 A 2 

4 


+ subleading 


where A is our UV regulator, m is our 1R regulator, A± is the area of the boundary, 0 is the 
noncommutativity parameter of the plane and a is the nonlocality scale of the dipole theory. 
The details of the expansion in ^ used to extract these leading divergences are discussed in 
Section 15.2.11 

In both cases, the contribution from these nonplanar diagrams does not have the right 
form to be interpreted as the sign of a volume law in the entanglement entropy and we must 
conclude that these nonlocal theories at least to first order in perturbation theory obey an 
area law. This can be contrasted with the strong coupling result which found clear signs of 
the volume law even for large regions [5]. Thus, the volume law must either only appear at 
higher orders in perturbation theory or it must require strong coupling. Consistent with our 
analysis, previous investigations of perturbative dynamics of the noncommutative theory p5j 
have shown that noncommutativity does not introduce any new perturbative UV divergences 
that cannot be reinterpreted as IR divergences. Thus, is it hard to see how the higher degree 
of divergence required for a volume law can arise in perturbation theory. We are lead to the 
conclusion that entanglement on distances of the nonlocality scale and volume laws require 
strong coupling and are not accessible to perturbation theory. 

The remainder of the paper is organised as follows: Section [2] describes the theories we 
study, Section [3] explains how the entanglement entropy can be computed perturbatively 
in these theories, Section [4] shows that the results for the free theory are unchanged in 
these nonlocal theories, Section 5T computes the first order correction in the coupling to 
the entanglement entropy in a real scalar 0 4 theory for a warm-up and for later reference. 
Section 5.2 extends the calculation to the real scalar on the noncommutative plane. Section 
5.3| reproduces the results for the previous two sections in the case of the charged scalar. 


Section 5A computes the result for the charged scalar in the dipole theory. Finally, Section 
|6] concludes with a discussion of these results. 


2 Theories 

The theories used in this paper are scalar field theories on M 1,3 where products of fields 
are replaced with a possibly noncommutative product denoted *. Three examples of this 
product will be used: the regular commutative one, the Moyal product associated with the 
noncommutative plane and the dipole product with a fixed nonlocality scale. See [TBJ for a 
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review of noncommutative field theory. The Euclidean action is 


Sr — d d x 


^<90 * <90(x) + ^ m 2 (j) * 0(x) + ^j-0 * 0 -k 0 -k 0(x 


The entanglement entropy in these three theories is calculated to leading order in the coupling 
A. The mass is present to serve as an IR regulator and will be taken to be small in the end. 

First, the standard commutative case, where (/ * g)(x) = f(x)g(x), is reviewed and 
presented in our notation in Sections [4] through 5T The entanglement entropy for this 


theory was studied in ra and the approach contained therein will be followed for each of 
the theories we consider. 


Second, in Section 5.2 the entanglement entropy of a held theory defined on the non¬ 


commutative plane, where 


(f*g)(x) = exp I -0^4-4- 
w F 1 2 <9^ <9(0 


f(x + £)g(x + t)\z=s= o, 


is studied. The noncommutativity is parametrised by the antisymmetric tensor 0. This 
theory has been studied perturbatively in pa]. In this case especially, the mass should be 
thought of as an IR regulator and taken to zero at the end of the calculation in order to see 
full effects of the UV/IR mixing present in this theory. We specialise to the case commonly 
referred to as the noncommutative plane where 0^" = 0 ( d lfJ, 6 2l/ — 5 2fJ, 5 lu ) for simplicity. 

Finally, the entanglement entropy of the a simpler nonlocal theory with a fixed nonlocality 
scale along a particular axis, known as a dipole theory, where 

(/ * g)(x) = f{x + ^ Q L {g))g{x - ^Q L (f)) 

is studied. 

This requires a global charge to identify with the dipole length. This can be obtained by 
using a charged scalar rather than the real scalar held theory we have discussed so far. The 
action for a complex scalar is 


Sr = / d d x 


-<90* * <90(x) + m 2 0^ * 4>(x) + ~-0^ * 0 * 0t * 0(x) + ^-0 t * 0 * 0 * 0 t (x) 


where there two 0 4 terms which are inequivalent due to our noncommutative product rad 
The result from the real scalar theory will be extended to this charged scalar theory in 
Section |5.3| , then the dipole theory will be studied in Section [54 


Setting Ql = aQu{ i), where a = ai is the fixed dipole length and i is a unit vector in the 
X\ direction, the terms in the action can be written in a more explicit form: 


J dx((fJ * 0)(x) = J dx(f)\x + -a)0(x + -a) = j dx(p\x)(f)(x), 

/ 111 1 

dx(jy{x + -a)0(x + -a)0(x — -a)0px — -a). 

Z/ Zj Zi Z 


( 1 ) 


4 These noncommutative products are constructed to ensure that integrals of products of fields are invari¬ 
ant under cyclic permutations. 
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In fact, renormalisability requires that we include terms X n f dx((j^4>)(x + \na)((jA <j))(x — 
|na) for all n id- However, the contributions from these terms can be obtained by simply 
substituting a —> na into the results for n — 1 and summing over n. The results in Section 
|5.4| are such that this sum is guaranteed to converge as long as the X n don’t grow too quickly. 
As the inclusion of these terms would not affect our conclusions, we will not consider them 
separately. 


3 Entanglement Entropy 

The standard technique of the replica trick is used to compute the entanglement entropy 
Id- This technique was used in a perturbative context in id. whose approach is followed 
here. 

Starting with p A , the reduced density matrix of the ground state of the theory in question 
for a region A, the idea is to evaluate 

S = -Tr(p A lnp A ) = — hiTr(/?^)| n =i, (2) 

by calculating Tr p n A for arbitrary n and analytically continuing. In this paper we will con¬ 
centrate on the simplest case where A is the half plane (A = {(aq, X 2 , x^) € M 3 |xi > 0}). 
The main result that will be needed can be lifted directly from mm- 

In Tr(p^) = In Z n — n In Z\, (3) 

where Z n is the partition function of the theory on an n-sheeted surface with a cut along 
the region A that connects the sheets. However, some details of this n-sheeted space will be 
needed in the argument to follow, so the rest of this section will define it more carefully. 


3.1 n-sheeted surfaces 

The density matrix can be written as a path integral, (at finite inverse temperature of /3) 

( fo \ p \< h ) = ( z ,)- 1 j V4,%T=tx~ s ‘ * 

where Z\ is a normalisation factor to ensure that Tr p = 1. Then the reduced density matrix 
for a region A is obtained by periodically identifying the held in the Euclidean time direction 
along A, the complement of A, while leaving the boundary condition along A untouched. To 
look at the ground state, /3 must be sent to infinity. We do this while keeping the cut along 
A near the origin. 

Then, 


Tr(p^) 


(Zi)~ n 



<t>(x£A,0 )=<p 2 


e 


—Se 



<t>(x£A ,0 )=4 >3 
(p(x£A,0+)=4>2 


e 


—Se 


J U( f ) 4>{xdA,0+)=<t> n 



This identification of boundary conditions can be replaced by defining the held theory on an 
n-sheeted surface with a cut along A that takes you from one sheet to the next. Calling this 
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n-sheeted surface (M. d \ A) n , the projection onto the sheet n : \ A ) n — s- \ A and the 

indicator function telling you if you are on the k th sheet Xk '■ \ A)” —> Zj, this means 

that <f> : (R d \ A) n -)■ R can be dehned as <F(x) = 0fc(7r(x))xfc(a0, so that 

I V$e~ SE , 

where Se for $ has the same form as that for each 0, since the action for each sheet is 
additive. 

With our simple region A, a half-plane, polar coordinates can be dehned in the x-r plane 
of \ A. Then the glueing required to create this n-sheeted surface is simply to identify 
6 = 2tt on one sheet to 9 = 0 on the next. Thus polar coordinates can be dehned on 
\ A ) n where 6 G [0,27rn), such that each interval of length 2n corresponds to a sheet, 
i.e. ir(r,6,y,z) = (r, 0 mod 2 t r,y,z) and Xk(r,0,V,z) = X[2n(k-i),2nk)(0)- 

This gives us the result from H3H3 cited above, as Z„ = f V&e Se . This path integral 
over $ is the path integral over the n-sheeted surface. 

4 Free Theory 

The hrst step is to understand the free theories where A = 0. The action for the free 
noncommutative and dipole theories is the same for that of the commutative theory, since 
the star product of 2 fields is the same as the regular product up to a total derivative floj . 
For the noncommutative theory, 

d„J(x) d vi ... d Un g(x) 

d/; 2 • • ■ d^fix) ... d Un g(x ) 

so that the quadratic term in the action is the same as for the commutative case up to a total 
derivative. As there are no boundaries, the only place this total derivative could make for 
a finite contribution is at the conical singularity introduced at the origin when considering 
the n-sheeted path integral. 

Around the origin this term contributes (note that the singularity is at the origin of the 
x-t plane and is not localised in the y-z directions), 

OO p 

lim A± / rd6 (-W Q'' 2 " 2 ... 0 , ' ni ' n <9 ... dn n (j) d Vl ... c0 n 0 ~ lim rd n (j)d n+1 q 6, 

r —>-0 ^ J / r —>0 ^ J 

n= 1 n 

where A± is the area of the y-z plane. As long as d n (f) <9 n+1 0 is regular at the origin this term 
will not contribute to the action. This means that 0 needs to be C°° at the origin, which is 
just the regular boundary condition imposed in the commutative case. 

For the dipole theory, direct calculation of the ^-product of two fields can be seen to 
reduce to the commutative result in Equation (JT|). 

Thus the free theory is the same for all three theories. 


d 4 x(f -k g)(x) = / d A x l - 0 ^ ... 0 ^" d, 


n =0 


'/U ■ ■ ■ 


= / gtx 


f(x)g(x) + d ril 


72=1 


Tr(pl) = {Z x )- n 
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4.1 Green’s functions 


Since the free theories are the same, they have the same Green’s functions. This Green’s 
function is straightforward in the polar coordinates introduced in Section 3T Since the 
action for $ living on the n-sheeted surface is the same as the action for </> living on any 
particular sheet, the local equation that the Green’s function must obey will be the same. 
The only difference is that 6 must be periodic with period 2i rn rather than the usual period 
of 2ir. The Green’s function for the held living on the n-sheeted surface is, from 


G n (x,x') = 


2nn J (27r) rfx 


d d ±p ± “ ro o 

2^ ak / 


k =0 


, Jk/n{.QT)Jk/n{.QJ ) n in nl\ / \ ivi (xi —x' ) / a\ 

dqq ' \ ' cos (k(9 - 0 ) n)e w±[x± x ±>, (4) 

q 2 + PI + rn 2 


where a 0 = 1, a^ 0 = 2, p± = (p y ,p z ) and x± = (x 2 , x 3 ). _L refers to the directions orthogonal 
to the cut introduced by the replica trick. 

The Euler-Maclaurin formula, 

x r /*00 "I 1 TD 

J2 a „F(k) = 2 / dkF(k) --F'(0)-2^ 7 -|i y Fra- 1 >(0), 
fc=o L>/o J 


can be applied to this Green’s function to replace the sum over k, 
G„M=r~ [ r d gq J ^ qr) f kiqr ] cos (h(0 - 


7T 

1 


(2ir) d± 

d d± p A 


127rn 2 J (2tt) c 


dqq 


E 

3> 1 


B■ 


2 j 


Jo 

d d± p A 


q 2 + p 2 + m 2 

dv[Jv{qr)J u {qr%= 0 ^ p± (x± _ x , ±) 

q 2 + p\ + m 2 

r °° {du) 2i ~ 1 [Ju{qr)Ju(qr') cos {v{0 - 9'))} u= o JPi{x[ _ x ^ ) 


(5) 


7rn 2j (2j)! J (27r) d± J 0 


dqq- 


T + P 1 


It will be useful to define G n {x,x'-,p) as 


G n (x,x'-, Py )=-L-J d £'£a„J 

J k =0 u 


rfgg i k,n ^ Jk, o n< " qr J 0 cos(fc (0 - e')/n)e iMx3 ~ x 3 )+^^-^) 


such that 


<T + Plj + P 2 + m 2 




( 6 ) 


G,.{x.:r’) = j -^G n (x,x'-,p s ) 

d d 

G n (x,x'-,p) = -- —G n (x,x'-,p) = ipG n (x,x'-,p). 


dx 


'2 LyU/ 2 
It is also useful to dehne f n (x,x') and f n (x,x';p) as 

f n (x,x') = G n {x,x') - Gi(x, x') 


(7) 


n 


2 -1 


127rn 2 J (27r) 
f n (x,x'-,p) = G n (x,x’-p) - Gi(x,x';p), 


££ r + 0 > i) 


g 2 + P]_ + m 


where G*i is the Green’s function on the 1-sheeted surface, that is just the regular Green’s 
function. 
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4.1.1 Single sheeted limit 


This Green’s function for the n-sheeted space must reduce to the regular Green’s function 
in the limit where n —> 1. Starting with our expression for the Green’s function in Equation 
Q, defining ip = 6 — 6' for convenience and setting n — 1, 


ri t 1 f dd± P ± 


E 

k =0 


[°° J Jk(qr)Jk{.qr') n ^ ip! ( x \ —x'i ) 
a k dqq— —--- cos{k<p)e p±{ ± . 

o q + p ± + m 2 


Equation (10.9.E2) in the DLMF [IE] provides a useful integral representation of the 
Bessel functions, which can be rewritten as, J n (z ) = f* p. e d zsm i~ n i) m Using this represen¬ 
tation and the fact that J- k (z ) = (—1 ) k J k {z), [^] 


OO OO /*7T 7 7 

J2a k Jk{qr)J k (qr')cos(k<p) = ^ / -^e i9(rsin7+r,sinK) - ifc(7+K) e^ 

k =0 k =-oo 


d'y 

2n ( 


t iq[r sin 7+r ; sin(</?—7)] 


Defining our position axes on the Xq-X\ plane such that x = (0, r) implies that x' = 
(— r' sin ip, r' cos <p ). Then defining q = (q cos 7, q sin 7), 


q- (x — x) 

OO 

^ a k Jk{qr)Jk{qr') cos (kip) 

k =0 


=q [r sin 7 + r' sin(<£> — 7)] 



iq-(x-x') 

2tt 


Finally, defining p = (q,p±), 


Gi(x, x') 


rfdp e ip(x-x r ) 
(27i) d p 2 + m 2 ’ 


which is the usual Euclidean Green’s function. 


4.2 Entanglement entropy in the free theory 

The entanglement entropy when A = 0 must be identical in the three theories as it was 
shown above that the quadratic terms in the action are the same. This can be seen more 
explicitly by using the approach from [TTj. Starting from Sa = — d n [In Z n — n In Z { \ n=v the 
part of the entanglement entropy which depends on the mass can be related to the Green’s 
function by 


d 

dm 2 


In Z n = 


1 

2 


d d x($ 2 (x)) n . 


5 Equation (10.4.El) in ,fB. 








In the commutative case, ( < h 2 (a:)) n = G n (x,x). In the non-commutative case, 


($*$(z)) n = exp 


i / d d d d 
2~ d( 2 d& dCi 


= exp 


dp y 


d d 


d d 


-0 

2 V di x d( 2 d£ 2 dCi 


mx + OHx + ())r 


G n {x + £, x + C) 


C=C=o 


2vr 


exp 


1 / 5 d 

2 p ’ V05i + SO 


?=C=o 

G n (x + £,X + C,Py) 


S=c=0 


dp v „ , 1„ . 1„ . N 

^G n (x+ -0p y i,x + -0p w i;p y ). 


That the ^-product turns out to just translate the argument of the Green’s function is an 
important theme of the calculation in this paper. 

The only difference for a complex scalar is that the mass term in the action is proportional 
to $1 * $ instead of <1> * $, however the expectation value of this leads to the same Green’s 
function and the same result follows. 

The dipole theory is identical except that translations by 0 times the momentum in the 
^-direction are replaced by translations by a. 

Thus, still for the non-commutative case, 


d 

dm 2 


In Z n = - 


1 

2 

1 

2 

1 

2 


d d x($ * <&(x)) n 



dp 

2tt 

dp 

2tt 


G n (x + -Opi, x + -Opi: p) 


G n (x,x;p) 



recovering explicitly the result from the commutative case by shifting the integration variable. 

However, this shift of the integration variable on the n-sheeted surface bears further 
investigation. It is sketched in Figure [l] 

This shift is well defined except for the region which gets translated into or out of the 
origin. However, this region has measure zero and cannot affect the result of the integral. As 
long as only a countable number of such shifts are done, these points can be omitted from the 
integral without changing the result. Finally, the integral over the whole n-sheeted surface 
can be written as a sum over the sheets and the Jacobian of this shift on each sheet is 1, so 
the Jacobian of the whole shift does not introduce any new factors into the integral. Thus 
shifting the variable of integration on this n-sheeted surface is allowed with no Jacobian, 
just as for the plane. 


5 First Order in Perturbation Theory 

5.1 Commutative theory 

We will start by computing the first order correction to the entanglement entropy for the 
commutative 0 4 theory. This was done previously in [T3], but will be repeated here with 
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Figure 1: Translations on each of the sheets of the n-sheeted surface (on the left) give a well 
defined map on the whole surface (shown for n — 2 in the polar coordinates described in 
Section 3.1 on the right), except for a measure zero set near the singularity at the origin. 


more explicit regulators that will allow a direct comparison to the nonlocal cases. From [H], 


In Z n = In j V(pe- SE[4,] 

— hi Z n fl — jy J d 4 x($ 4 (x)) 0 + ... 

= In Z nfi - ^ J^d 4 x [G n (x,x )] 2 + ..., ( 8 ) 

where J n denotes integration over the n-sheeted surface and In Z n ^ is the k th order term in 
a A expansion of In Z n . Generally, adding subscript will denote the order of a term in a A 
expansion, e.g. X = A" 0 + X x + X 2 + ... 

The entanglement entropy can be calculated using Equations (j2]) and (|3]), 


In Tr {pa)i= In Z n , x - n In Z x ,\ 

= - J d 4 x [G n (x,x)] 2 + J d 4 x [Gi(ic,ic)] 2 

q\ r 

= - — ^ d A X [2Gi(x, x)f n (x, x) + f 2 (x,x)\ . 


Recalling from Equation (JtJ) , 


fn(x,x') = 


n 2 — 1 f d d± p J 

127 m 2 J (2n) d 


dqq 


d v [Ju(qr)Ju(qr')\u =o 
q 2 + p 2 ± + m 2 


e ipj.(*x-*i) + (j > i). 


The j > 1 terms don’t contribute m so they will be dropped in what follows. This is the 
same on each sheet, so the integral over the n-sheeted surface is n times in integral on one 
sheet. Finally, fi(x,x') = 0, so d n f 2 (x,x')\ n= i = 0 and 


fiA f 

Si = -d n [InTr (^) 1 ] n=1 = yjj / d 4 xG x (x, x)d n [nf n (x, x)] n=1 


(9) 
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Sl = ELL. f riTd4 , f d ypy*P‘ 1 r dm MMsrWf })--» . 

12n ■ 4! J J (27r) 6 k 2 + m 2 J 0 q 2 +p 2 +p 2 + m 2 

Schwinger parameters are introduced to allow the denominators to be combined, using 

dae~ Aa . 

1 

This allows us to regulate the UV divergence in S\ by introducing a factor of e «a 2 . as 
was done in previous perturbative studies of noncommutative theories [T5]. This regulator is 
convenient in the noncommutative case and is used here so that the results can be compared. 
Using Equation (25) from p.146 in volume I of | [TUj . 



dte- pt ~% = J-KAJdp), 
V P 


the effect of this regulator is 


dae ~“ p2 -^ - 2 


31 = i-K 1 X 

A p \ A 


2 __ 2p 

- f A 


A p- 


£->Q 1 

p 2 ' 


( 10 ) 


Thus it regulates the UV and leaves the IR unaffected. This can be seen simply from the 

fact that e~~^ vanishes for a <C A -2 and goes to one for a A~ 2 . A mass m regulates the 

2 

IR by contributing a factor of e~ am , which has the opposite behaviour. 

Introducing these Schwinger parameters and regulating, 


Si = 


XA± 
3 • 2 3 


dr d*kdp dp z dq 


(2vr 


All the momenta integrals except q are Gaussian, 

Si = g A ^ r 3 J drdq J dad/3^e~ Pq ~ am ~^~ Pm ~~^d u [J u (qr)J u (qr)} u=0 . 

Using Equation (10.22.E67) from the Digital Library of Mathematical Functions (DLMF) 

02 ], 


oo 1 (a 2 +b 2 ) / nh 

te- p t Uat)J v {bt)dt = , 


( 11 ) 


the q integral can be evaluated. This along with the fact that d u I u (z)\ u= o = —Kq{z^ gives 

XA i ' '■°° - ~ 2 - ’ - - 1 / - 2 


Si = - 


3 ■ 2 10 tt 3 


dr / dadfd 


r — XX* — rvm? -g*-— Rm? -L- - - I T“ \ 


a 2 ft 2 


2 
e 


^ A ° 2/3 ) 


0 Equation (10.38.E4) in the DLMF 
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Equation (21) on p. 131 of [19], 


arccos(-) 


vV - « 2 

after substituting r 2 —* t and setting a = y = gives 


dte at K 0 (ty ) 


f-n 1 
J! —> 

y 


Si 


XA± 

3 • 2 10 tt 3 





a' 


e “ m 



P 


Looking at the a integral first, 




—2mAK\ 


mP 1 _ a. 
a A? 


7?-°. A 2 


( 12 ) 


by substituting a —^ in the first line and using Equation (10) as well as in the second. 
This recovers the A 2 divergence seen previously in this case m 
Using Equation (29) from Volume 1, p. 146 of [T9] 


the f3 integral gives, 


t »-i e -rt-itdt =2 ( — ) K v {y/ap) 


df3 -p(2m\ 01 2m A 

P \ A J A fan 2 


as Kq(z) —> — In z as z —> 0. This reproduced the logarithmic divergence seen previously in 
this case [14j and makes explicit its form in our regularisation scheme. 

Combining, the first order in A correction to the entanglement entropy in the commutative 
theory is 


Si, 


Comm. 


= — 3A 


A i A 2 


In 


A 2 


3 2 • 2 10 7t 3 4m 2 


(13) 


This is proportional to the area of the boundary of A, that is A±, and the leading divergence 
is of order A 2 , so this result fits with the area law picture discussed in the introduction. 


5.2 Noncommutative theory 

Next we will compute the first order correction to the entanglement entropy for the noncom¬ 
mutative 0 4 theory. Similarly to the commutative theory, 


In Z n = In / T>cf)e 


—Sb[4>] 


— hi Z n fo — /’d 4 z($*$*$*$(: C )}o + .. 
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Using the associativity of the ^-product, this can be written as 


d 4 x(& * * <f> * $(a;)) 0 = / d 4 x ( exp 


^exp 

>( 

^exp 

X 

^exp 

X 


d 

d 

d d \ 

% 

1 

d&dCi) 

d 

d 

d d \ 

drji 

0^2 

drj 2 d$i) 

d 

d 

d d N 

dli 

dn 2 

di2 dhi\ , 


?=C=o 


•n=s =o 


7=k=0 


{®(x + £ + rj)$(x + £ + + c + l)$(x + c + «))■ 

The usual Wick’s Theorem can be applied to calculate the four-point function, 

($(w)$(x)$(y)$(z)) = G n (w, x)G n (y, z) + G n (w,y)G n (x, z) + G n (w, z)G n (x,y). (14) 

The key point is that while the conical singularity breaks the translational invariance 
in the xq-x\ plane, it is preserved in the ^-direction. Thus the star product reduces to 
a translation in the Xi-direction by an amount determined by the momentum in the re¬ 
direction. Defining G n (w,z) = f ^-G n (w, z;p y ) as in Equation (|6]), 


? <9 d 

exp 1 2 e d^ar 2 


dp f 1 d 

G n (w, z) — I -^exp ^-p y G-Q—)G n (w,z\p y ) 


= I ^G n ( W + ±p y Oi,z;p y ), 


this can be used to evaluate the 4-point function, 

/ d 4 x < * <h * $(rc) >o 

J n 

= [ d 4 x 


dk y dp. 


irL ' v G n (x + ^Qk y i,x + ^Qkyi] k y )G n (x + \>Qp y i,x + ^®p y v,p y ) 


(2n] 

+ G n (x + -Qkyi, x + -®(ky + 2 p y )% ky)G n (x + -0(2 k y + p y )i, x + -Qp y r, p y ) 
+ G n (x + ^Qkyi, x + ^®k y i ; k y )G n (x + ^0(2 k y + p y )i, x + ^0(2 k y + p y )i;p y ) 


Then, by shifting the spatial integral, 
r r dkydpy 


= / d 4 x 


(27t) 5 


G n (x, x ; k y )G n (x + ^Q{p y - k y )i, x + ^Q{p y - k y )r,p y ) 


+ G n {x - ^®p y h x + ^®Pyt; k y )G n {x + ^ ®k y i , x - ^®k y v, p y ) 

+ G n (rc, rc; k y )G n (x + ^ ®(k y + p y )i, x + ^ ®(k y + p y )i\ p y ) 

In [15] it is seen that the effects of the non-commutativity manifest themselves in the 
diagrams where lines cross each other. This is also present here, as Figure [2] shows that it is 
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Figure 2: Vacuum bubble diagrams at leading order in a real scalar A0 4 theory. The only 
vacuum bubble where lines cross is the second one. This is the only one which is affected by 
the non-commutativity, as discussed in [ T5] . 


only the second term that involves lines crossing. The other two terms are two self-coincident 
Green’s functions - the same result as was found in the commutative case in Section 15. II and 
pH ], The second term, which corresponds to the nonplanar diagram, is the only one which 
is different than what was found in the commutative case. 

The entanglement entropy can be calculated using Equation (J2]), 


Si = -d n [In Z nA - n In Z i,i ] n=1 

= Tf®”( / d“x 2Gi(x,x)n/„(x,x) + 


(15) 


dkydPy 

(2tt) 2 


G i 0 , x + Qpyi] k y )nf n (x, x - Qk y i; p y ) 


where the fact that the spatial integral can be shifted, that the momenta can be renamed, 
that Gi(x,x]p y ) = G\(x + a,x + a\p y ), that f n (x,x' ,p y ) = f n (x,x — p y ) as long as x 2 = x ' 2 
and that /i = 0 so that the terms with f 2 can be ignored have all been used. The j > 1 
terms in f n have also been dropped again, which allows us here to write the integral over the 
n-sheeted surface as n times the integral over a sheet. In the commutative case, it was clear 
that these j > 1 terms do not contribute [12] . In Appendix [A] it is argued that the leading 
divergence must be entirely contained in the j = 1 term even in this noncommutative theory. 


5.2.1 New contribution from the nonplanar diagram 


The first term in Equation (15) is the contribution from the two planar diagrams. These 
give the same result as in the commutative case, namely 21032^3 l n from each diagram. 
However, the nonplanar diagram gives a new contribution to the entanglement entropy from 
the non-commutativity. The contribution from this nonplanar diagram will be denoted 

“Snonplanar, 



Gi(x, x + Gpyi: k y )d n [nf n (x, x - Qk y kp y )] n=1 

d A kdp y dp z e iekxp y f°° d u [J u {qr)J u (qr ')} v=0 

(27t) 6 k 2 + m‘ 2 J 0 ^ q 2 + p y + pi + m 2 


1 


where r ' 2 = (r*— Qk y i ) 2 = r 2 + {Qk y ) 2 — 2Qrk y cos(p and A± is the area of the X 2 -X 3 plane 
that bounds the region for which the entanglement entropy is being calculated. 

The next step is to introduce Schwinger parameters and to regulate this integral in the 
same manner as the integrals for other perturbative calculations in this noncommutative 
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theory were regulated in [13] . as discussed in Section 5.1 

XA i 


5, 


nonplanar 


2 3 3% 


drdcj) 


d A kdp y dp z 
(27r) 6 


dq / dad^qTe- ak2 < q2+p2 ^^ am2 -^ m2 


e iekxPy d u [J u (qr)J u (qr')\ u=0 . 

The p y , p z and k except for k y integrals are all Gaussian (recall that r' is a function of 


v yJ-> 


S, 


XA 


nonplanar 


2 8 3 2 7T2 


drdcpdkydq / dad/3 


qr 


o ay/^yjA a(3 + @ 2 


m ~P m ' 


^ d u [Ju{qr)Ju(qr')) v =o- 

In order to make explicit some of the symmetry between r and r', p and p can be dehned 
such that r = p sin p and k y = ^cosp, with p G [0, oo) and p G [0, 7r] . Then defining 
g((j), p) = v^^Fsm^^cos^, gives r’ = pg((j), p) in these variables. Performing this change of 
variables, 


5, 


XA i 


nonplanar 


2 8 3 2 7rf0 


d u \ v= Q / dpdpd(f)dqdadj3 


qp sm p 


Aa/3 + 0 2 




From the DLMF (10.22.E67) [IS], 


9 9 1 (a 2 +b 2 ) / nh 

te-^UatiUbt) = —e— &~I U ( — 


(16) 


so that, 


5, 


XA i 


nonplanar 


2 9 3 2 tt20 


d„\ v=0 / dpdpd<pdadfd 


p sm p 


a/3z y/4 a/3 + @ 2 


_^p2 cos 2 y- p 2 sm 2 V +g\ V ,<t>) - l_ am 2_l_p m 2 ( p 2 

e 4,3 H I v l—g[(j),p)smp 

Now p and a can be rescaled to simplify this expression as p —* 2y/~/3p and a —> 


S, 


XAj 


nonplanar 


2 6 3%i0 2 


d v \ v=G / dpdpdtfida 


p sm p 


a 


\J ot + 1 


dj3e 


-p(r. 


l2 +e4^ 




e -a/co.V-p , [.inV+9 1 («,*■)]/__ ( 2 / 9 (^ if) sin p) . 

Equation (25) from p.146 in volume I of [T9] , 


dte- pt ~% = J-KAy/fip), 
V p 


(17) 
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allows the /3 integral to be evaluated, 


dj3e 


-p(r 


"e 2 A 2 , 


iBO 


4 7 — 


^ + 0 2 m 2 a 

™ 2 + e4ui 

/ 4 


Ah 


— + @ 2 m 2 ct ] ( m 2 + 




\0A 2 a/« 


A 2 

+ Om 2 \/a ) . 


)(• 


4 ) 

0 2 A 2 a 7 


Using the identity (9, y | y=0 /y(^) = — A' 0 (^), 


A, 


AA 


nonplanar 


2 6 3 2 7r 2 0 


dpdpd^da f s my; e~ p2 [“ cos2 ¥,+sin2 ^ +s2 

yCTy Ct + 1 




AT 0 (2p 2 g((j), p) sin </?) A', + @m2 


Taking a large A limit of this expression and expanding K\ (x) « A for x —)■ 0 allows 
us to extract an overall quadratic divergence. However, more progress can still be made by 
evaluating the p integral. 

Using in order Equation (23) from p. 131 of [19] and (15.9.E19) of [ IS] . 


dpp 2 e Ap ~K 0 (Bp 2 ) 



dx\fxe Ax K 0 (Bx ) 


2 r( 2 )(A + A)f 

3 1 

7T 2 1 


8 V 2 B 


if -1 


T 2 A 1 ( -, 2; 




A- HA 

a + aJ 


where P\(x) is 
2 

degree. 

Dehning z = 


the appropriate branch of the associated Legendre function with non-integer 
aC ° s2 and recallin § that p) = \/l + sin 2 p cos 0, 


AA ± r , G(a) „ ( 4 

nonplanar 31132^0 J Q ^ ^^+1 ' \QA^ 

rn r- 2 -K 1 P\ (z) 

G(a ) — dip d(j) - 3 - 2 , 

Jo Jo [g((p, p)]^ y/^nip V z 2 — 1 



and 


where G(a) is dimensionless and finite for a G (0, oo). 

At this point, the asymptotic behaviour of G(a) can be analysed numerically, as no 
analytic formula for this integral was found in the tables consulted. However, while analysing 
this asymptotic behaviour, we found that G(a) = gives an exact match up to high 

numerical accuracy across the many orders of magnitude that were checked]^] 

'The only potential divergences in the integral for Nonplanar come from the regions of small and large 
a. If the reader is uncomfortable with this numeric argument, this functional form for G{a ) could also be 
thought of more conservatively as a function with the right asymptotic behaviour to reproduce the correct 
divergences in this integral. 
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Using this result for G(a), 


S, 


XA 


da 


nonplanar 


2 7 3% 3 0 


-K 


A, 


a a + 1 \ 0 A 2 A/a 


+ Qm 2 ^/a 


Note that this result is invariant under @A 2 -h- 0m 2 , another sign of the UV/IR connection 
in non-commutative theories. 

This integral has two regulators, A and m. The only other dimensionful parameter is 0, 
so the only dimensionless products of these regulators are ^ and 0mA. As is familiar from 
the UV/IR mixing in this theory, the limits A —» oo and m —> 0 do not commute. This can 
be resolved by taking ^ —> 0 while fixing 0mA. Then taking the limit m —> 0 or A —» oo 
first corresponds to the limits 0mA —» 0 or 0mA —» oo respectively. [^] 

Introducing 7 = \fa, 


5, 


XAi 


nonplanar 




2 6 3 2 7t 3 0 J 0 
XA ± A 

2 8 3 2 7r 3 0m J 
AA, A 2 - 1,1 


dy 


7 2 + 1 


*■' ¥ 


QmAy 


+ 


QmAy 


dj 


7 2 + 1 


GmA'y 


+ 


GmA'y 


Q 2 m 2 A 2 \ 

4 J 


2932^3 _ 9 2 m 2 A 2 


^_L _ l 11 ( 0 2 m 2 A 2 ) 

2 7 3 2 7r 3 0 2 m 2 1 


0 2 m 2 A 2 


(18) 


where the last line uses Equation (2) from Volume 2 p.216 of [15]. 

This result illustrates the UV/IR connection in non-commutative theories. If the IR 
regulator is removed first (0mA <C 1), Skonpianar ~ Aj_A 2 - a quadratic UV divergence. 
However if the UV regulator is removed first (0mA 1), Skonpianar ~ qAx ’ allowing the 
same divergence to be interpreted as an IR divergence. In addition, whether e ~™ 2A ~ is taken 
to be large or small there is a logarithmic divergence as is found in the commutative case. 
However, here there is the additional option of keeping both regulators, that is keeping 
(0mA finite, which eliminates the logarithmic divergence seen in the commutative casej^ I 11 
particular, there is a natural choice of IR regulatoiJ^J m = Ah where 

C( XA ± XA ± A 2 

nonplanar 2 7 3 2 7 T 3 0 2 m 2 2 9 3% 3 ' 


From a mathematical point of view, this UV/IR connection can be seen to originate 
from the translation of the arguments of the Green’s function. In the commutative theory, 
Snonpianar ~ f n dxG n (x,x)f n (x, x) where as in the noncommutative theory, the non-planar 
diagram made a contribution of the form SAonpianar ~ J dxG n (x, x + Qp)f n (x, x + Qp). If an 
IR regulator is imposed, this momentum cannot vanish and regulates the integral. This can 
be seen more clearly in the dipole theory (analysed in Section 5.4) where the fixed translation 
regulates the UV divergence of the integral. 


8 This discussion applies even if we want to think of to as a physical mass, as the ratio ^ will still vanish 
if to is fixed while A —> 00 . This case corresponds to ©toA —> 00 . 

9 Note that if a 0 —> 0 limit is taken, this option is no longer available and the commutative result is 
recovered, although the exact form of the logarithmic divergence depends on how the 0 limit is taken. 

10 See Section 6 of [TS] 
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Figure 3: Vacuum bubble diagrams at leading order in the noncommutative charged scalar 
A0 4 theory. The two on the left come from the Ao0' * <fi * (fy * cj) term in the action whereas 
the two on the right from the Ai $ * cj) -k 0 * term. 


It is important to note that contributions from the j > 1 terms in Equation ([5]) were 
dropped at the start of this section and are not present in Equation (18) or elsewhere in these 
results. However, as is discussed in Appendix [Aj these do not affect the leading divergence 
in SAonpianar or the conclusion that there is no volume law. 

In contrast to strong coupling results, which saw signs of a volume law for the entan¬ 
glement entropy even with large regions, this perturbative calculation is only sees an area 
law. The leading divergence in SAonpianar is quadratic and proportional to the area of the 
boundary of the region, A±, in line with the area law discussed in the introduction. 


5.3 Charged Scalar 

For the commutative and the noncommutative theory, the difference when considering a 


charged scalar comes in at the Wick contraction in Equations 
the real scalar 


and (14) respectively. For 


A {(p(w)(j)(x)(j)(y)(j)(z)) = A (G n (w, x)G n (y, z) + G n (w, y)G n (x, z ) + G n (w, z)G n (x, y )), 
whereas for the charged scalar this must be replaced with 
Ao {w)<t>(x)<fi (y)<j>(z)) + Ai (<f>\w)<j>{x)<l>{y)<ft{z)) 

= A 0 {G n (w, x)G n (y, z) + G n (w, z)G n (x, y)) + Ai (G n (w, x)G n {z, y) + G n (w, y)G n (z, x )). 

In the commutative theory, the helds in the 4-point function are all inserted at the same 
point, that is w — x — y — z. Taking into account the difference in the normalisation of the 
(j) 4 term in the action, the only change is to replace an overall factor of Ap by xhis 

has no effect on the intermediate steps of the calculation and can just be carried through 
straight to the final result: 


Si 


,Comm. ^ 


(Ao + Ai)Aj_A 2 A 2 
3 • 2% 3 4m 2 


For the noncommutative theory, it is a simple matter of writing out the ^-products 
explicitly and following through similar transformations of the integration variables as in the 
previous section. This procedure gives 2Ao + Ai times the commutative result plus Ai times 
the result for the nonplanar diagram already encountered for the real scalar. This result can 
be obtained directly by looking at the 4 diagrams in Figure [3] and realising that the only 
term proportional to Ai can give a nonplanar diagram. 
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Thus the result for the noncommutative theory with a charged scalar is 


A A 2 A 2 

<Si,nc —t(2A 0 + Ai) -—In ■ 


3 ■ 2 9 7t 3 4m 2 


-A, 


A ± A 2 
3 ■ 2 8 vr 3 


-lu^eh^j 

-| Q 2 m 2 A 2 

J- 4 


5.4 Dipole theory 

For the dipole theory, the explicit form of the interaction terms was written out in Equation 
Q. Thus, 

In Z n = In J V(pe- SE[4l] 

= In Z nfi — f d 4 x /^ < h t (a;) < f ) (a;) < f )t (a;) < f)(a;) + + ^-a)$(x + — ^a)<I>i(;r — ^a)\ + 

Jn \ 4 4 2 2 2 2/ 0 


Applying Wick’s Theorem, using the facts that Gi(x, x) = G\(x + a, x + a) and f n (x + 
a, x ) = f n (x, x + a) (when ignoring the j > 1 terms) and shifting the integral, 


In Z Ut i — — - / d 4 x[2A 0 G n (x,x)G n (x,x) 

^ J n 

+ Ai ^ 'G n (x +^a,x+ ^a)G n (x - ^a,x - ^a) + G n (x + ^a, x - i a)G n (x - ^a, x + ]xi) 
Si=-d n [In Z n> 1 - n In ^i,i] n=1 


= -d n 


d 4 x 


(2A 0 + Ai)G*i(a:, x)nf n (x, x) + AiGi(a;, x + a)nf n (x, x — a) 


n=l 


Again this is as expected from the diagrammatic approach. Only the single nonplanar 
diagram gives a new contribution and the 3 planar diagrams give contributions identical to 
those in the commutative theory. 

Focusing on the contribution from the nonplanar diagram, the explicit forms of G\ and 
fn give 


5, 


nonplanar 


4AAj_ 

12vr-4 


rdrdcf) 


d A kdp v dp z e lkxa 


(27t) 6 k 2 + m 2 J o 


' dv[Jv{qr)Ju{qr’)\y =0 

99 q 2 +p 2 v +p 2 z + m 2 ’ 


where now r /2 = (f— ai) 2 = r 2 + a 2 — 2racos0. 

Introducing Schwinger parameters and regulating, 


S, 


nonplanar 


AA ± 

2 2 3tt 


drd(J) 


d 4 kdp y dp z 
(27r) 6 


dq 


dadf3qre~ ak2 -P [ q2+p y +p2z ] ~ ^~ am2 ~jh ~ /3m2 
e lkxa d v [J u (qr)J u (qr’) } u=0 . 


In this case, all the momenta integrals except q are Gaussian, 


5, 


AA, 


nonplanar 


2 8 3tt 4 


drdcj)dqdadf3 


q r — f-G — am 2 —0q 2 — —hr — 0 m 

__ p 4a cA 2 ^ /3A 2 ^ 


a 2 /3 


^ m ^[^(gr)J I/ (g7- , )] v=0 . 
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The a integral can be factored out to give, using Equation (10), 




2 m 



Ki 


A—Kxi 4 ttt , , 

-y —Ki (ma) 

a 

m— >0 4 

-> — 


or 



This factor came from evaluating £0(0, a£) which goes as ~ \ as expected. The fixed 
nonlocality scale has regulated the UV divergence in this case. In the dipole theory the 
distance of the translation is fixed, as opposed to the non-commutative case where the 
translation is proportional to the momentum in the ^/-direction which can vanish in the IR. 

Using Equation ([TTj), 


5, 


nonplanar 


XA ± 
2 7 37r 4 a 2 





rr 2 , /2 

I _ r ~f~ v 

d<t>p2 e 4/5 K o 



e ^ 2 


2 


Rescaling r —y ar and /3 —y a 2 f3 to make them dimensionless (r / —y a\Jr 2 + 1 — 2r cos < 
under this) and defining H(/3) as the part of the previous equation enclosed in brackets, 


5, 


XA\ 


nonplanar 


2 7 3vr 4 a 2 


df3H^)e~^^ a2m2 . 


H(f3) is dimensionless and finite for /3 e (0,oo). The integrand is exponentially sup¬ 
pressed for small (5 and numerical evaluation of the r and 0 integrals confirm that H(/3) ~^°> 

0. The other potential source of a divergence is at large /3 and numerical integration finds 
that H(/3) / ^°°> y leading to a logarithmic divergence at large /3 that must be regulated 
by e ~P°- 2m \ 



ln(a 2 m 2 ) + 0 (m°), 


to leading order in the small m limit. 

Thus S'nonpianar has only an IR divergence in the dipole theory. The leading divergence in 
the j — 1 term is 


5, 


XAj 


nonplanar 


3 • 2 6 7r 3 a 2 


[— ln(a 2 m 2 )] , 


(19) 


however there will be contributions to this order from the j > 1 terms which were dropped. 
The the conclusion of this analysis is that the nonplanar diagram does not contribute to the 
leading divergence of entanglement entropy at this order as it is subleading to the contribu¬ 
tion from the planar diagram. 

The nonlocality introduced in the dipole theory does not affect the area law, as the 
total entanglement entropy at this order in perturbation theory is dominated by the planar 


20 



















diagrams which matched the result from the commutative theory. Even the subleading 
terms we have analysed do not follow any sort of volume law as they are not proportional to 
the lcngthscale of the nonlocality. The only effect of the nonlocality is to regulate the UV 
divergence otherwise present. Similar behaviour was observed in ra, where one of the ways 
that the nonlocality manifested itself was by softening divergences in nonplanar diagrams. 

6 Final remarks 

In this paper we computed the first perturbative correction to the entanglement entropy in 
two nonlocal theories, a 0 4 theory defined on the noncommutative plane and a dipole theory. 

The contribution to the entanglement entropy in each of these theories at first order in 
coupling comes from vacuum bubble diagrams. The planar diagrams give the same contri¬ 
bution in all three theories. However, the nonplanar diagram is affected by the modified 
^-product. Never the less, these diagrams do not modify the area observed in the commuta¬ 
tive theory. Thus, at this order in perturbation theory and for the region considered at least, 
all these theories follow an area law with no sign of a volume law, as opposed to the strongly 
coupled case where the signature of the volume law could be seen even for large regions. 

In the commutative theory it has been shown that the modification to the entanglement 
entropy at first order in perturbation theory can be absorbed into the renormalisation of the 
mass pa- It would be interesting to see if a similar interpretation can be made in the case 
of the theories considered here. 

Finally, a comment about the commutative limit. Since the quantities dealt with in the 
paper are not UV finite, this is not a straightforward issue. The general pattern is that the 
nonlocality has served as an additional regulator that softens certain divergences. Thus, if 
the nonlocality is removed, these divergences reappear and the commutative limit applied 
to the final results is not smooth. 
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A Analysis of the potential divergences from the j > 1 
terms 

This analysis follows that of ra. where it is found that the leading divergence when the 
Green’s function is evaluated at coincident points is entirely contained in the j = 1 term. 
The Green’s function for the scalar field on the n-sheeted space was given in Equation 
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G n (x,x') = — 


dk 

71 

1 


d d± p ± , Jk(qr)Jk(qr') n ta iVl (x . _ x > 

77T^7 / d( M 0,2 1 2 cos ( 6 ~ &))e ip±{x± x±) 

{ 2 n) d± J 0 q-+p z ± + m 2 
f d d ^p L 


127rn 2 J (2n) c 

fYi ™ 2 j ( 2 J) ] j 


Jo 

d d± p A 

(2n) d - 


' dqq dAMqr)Mqr%=0 c i P I (x, -»/ ) 

g 2 + + m 2 


dqq ^ 2i 1 ^ qr ^ J ^ qrn> cos ( u ( 9 ~ 9 '))]y=o c » Pl (x, ) 

g 2 + + m 2 


The first term is independent of rt and did not enter into the calculation of the entanglement 
entropy. The second term was the subject of our investigation. However, the third term was 
dropped with the claim that it could not introduce any new divergences. This appendix will 
justify this claim. 

We start by revisiting the entanglement entropy in the commutative theory. In this case 
from Equation (J9]) 


S ~ 


J rdrGi(r , r)f n (r,r), 


where only the contributions to the divergences in the final result have been kept. 
The Green’s function when evaluated at coincident points gives a A 2 divergence 


Gi(r,r) ~ / d 4 p 


j p 2 + m 2 

~A 2 — m 2 log A 2 . 
The f n term has the form 

du[Ju(qr)J u (qr)\u =o 


fn(r,r ) ~ / d 2 p ± qdq 


q 2 + p\ + m 2 


da p 3 dp e 


— a(p 2 +m 2 ) — - 


Y / d 2 p ± qdq 
j> i J 


l lJv(q r )Jy{.qr)\v=o 
q 2 + p 2 + m 2 


The momentum integrals can be evaluated when the function is evaluated at coincident 
points 


d 2 p L qdq 


~ e 


- ir Jjlr 2 


Jy{qr)J v (gr) 
q 2 + p\ + rn 2 

log A 2 . 


= J d/3 pdp qdq J u (qr)J u (qr)e +p +m > 


1 



This must be integrated over r 



rdre 


_i r 2 

2 h, 



1 

V 2 e 


v + 0(e), 


where a small e has been added to regulate the integral. It is only divergent because dj* 1 
was passed though the integral sign. Once this derivative is applied, e can be safely taken 
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to zero. A calculation of terms 0(A°) would require a more careful analysis, but this is 
sufficient for extracting the leading 0(logA 2 ) divergence. Thus 


[ d 4 xf n (x,x ) ~ A± log A 2 


= A± log A 2 

-i + £° 

/ 

j > 1 


j > 1 - 


This shows that all the j > 1 terms vanish when the Green’s function is evaluated at 
coincident points and the divergence is entirely contained in the j = 1 term. 

In the noncommutative and dipole theories, the Green’s functions are evaluated at points 
separated by the length scale of the nonlocality rather than at coincident points. However, 
we saw that the source of divergences in the entanglement entropy was regions in the integral 
where this separation vanishes. This analysis shows that these divergences are contained in 
the j = 1 term. 
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